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A parameterization of the density operator, a coherence vector representation, which uses a basis 
of orthogonal, traceless, Ifermitian matrices is discussed. Using this parameterization we find the 
region of permissible vectors which represent a density operator. The inequalities which specify the 
region are shown to involve the Casimir invariants of the group. In particular cases, this allows the 
. . determination of degeneracies in the spectrum of the operator. The identification of the Casimir 

CO . invariants also provides a method of constructing quantities which are invariant under local unitary 

' operations. Several examples are given which illustrate the constraints provided by the positivity 

' requirements and the utility of the coherence vector parameterization. 
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I. INTRODUCTION 



" ■ The density operator must satisfy three important requirements. 1) It must be Hermitian. 2) The trace of the 
_ -, , density operator, when appropriately normalized, must be one. 3) It must be positive semi-definite. The third of 
^ ' these requirements has been found to be vital in quantum information theory, and in quantum mechanics itself 

[ Perhaps the most important place this has arisen is in the identification of positive and completely positive maps 
p\| which can be used to identify entangled states 0)0 and to classify quantum channels (0,0 and references therein). 

For both of these problems, but in particular the latter, a parameterization of the density operator is often useful. 
04 : This provides an explicit way in which to identify when the channel is unital, trace preserving, and/or completely 
' positive (see for example 0). In addition, positivity requirements place restrictions on physically realizable quantum 
[ transformations 

■»«^ Here we represent the density operator using a basis of orthogonal, traceless, Hermitian matrices. This reprcscn- 

i-Q tation is the generalization of the Bloch or Coherence vector for two-state systems which is commonly used (see 8] ) . 
Oh! While the geometry of the space of density operators for two-state systems is relatively simple, the geometry of the 
' space of density operators for higher dimensional systems is considerably more complicated. The positivity (or more 
■ precisely, positive semi-definiteness) conditions are therefore more difficult to express succinctly for higher dimensional 
^ ' systems. The inequalities given in this paper give necessary and sufficient conditions for a Hermitian operator to be 
positive semidefinite. 

IL^ This set of inequalities can be expressed in terms of a distinguished set of unitary invariants, the Casimir invariants. 

This is a particularly notable relationship since the Casimir invariants are associated with the "good" quantum 
numbers of a quantum system and thus have direct physical interpretation. They specify the set of quantities 
5— i which are invariant under a given set of unitary transformations. This has found many important applications for 
9^ I modelling of physical systems, and more recently, in quantum control of spin systems [10| . In addition, the Casimir 
invariants and positivity requirements are expressed in terms of the coefficients of the characteristic polynomials. 
These coefficients, and their ratios, were found to be entanglement monotones Entanglement monotones could 
provide some insight into the problem of finding suitable entanglement measures since they satisfy an important 
requirement of such measures; they do not increase, on average, under local operations and classical communication 

This paper can be divided into three main parts (excluding the Introduction and Conclusion). The first part gives 
the generalized coherence vector representation of the density operator and the Casimir invariants in terms of the 
coherence vector. The second part gives positivity conditions for the density operator in terms of the trace invariants 
as well as the coherence vector. The third part gives some examples of the utility of the structures presented in the 
first two parts. 



II. COHERENCE VECTOR/CASIMIR INVARIANTS 

In this section we present a coherence vector representation for an A^-state system with particular normalization 
relationships which differ, for example, from 8}- This is the generalization of the Bloch sphere representation for two- 
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state systems. The coherence vector, in our parameterization, has unit magnitude for pure states and has magnitude 
strictly less than one for mixed states. We will then show how to construct the Casimir invariants of the system in 
this parameterization. Using a completely analogous construction, we are able to provide a distinguished set of local 
unitary invariants for composite quantum systems. 



A. Pure States in A'^-Dimensions 



Any density operator can be expanded in any basis of orthogonal, traceless, Hermitian matrices. Here we adhere 
to the following conventions. We will use the following normalization condition for the elements of the Lie algebra of 
SU{N) 

Tr(A,A,) = 2%. (1) 
We will also choose the following relations for commutation and anticommutation relations: 

[Xi, Xj] = 2ifijkXk (2) 

and 

4 

{Xi,Xj} ^ —Sijl + 2dijkXk, (3) 

where the fijk are the structure constants and the dijk are the components of the totally symmetric "d— tensor." 
These two equations may be combined more succinctly as 
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XiXj = -^^ij + ifijkXk + dijkXk- (4) 



Using these conventions, we may express a pure state for an N x N density operator as 

(5) 




This representation is called a coherence vector representation with n the coherence vector. The constant is a 
convenient one such that for pure states 

n-n=l, and n-kn = n, (6) 

where the "star" product is defined by 



r^ NiN-1) 1 , , 
(a * b)k = y j^—^ diJka^bj. (7) 

This can be proved by direct computation using Eq. 

Orthogonal pure states, e.g., |ai) and |a2) with corresponding density operators pi — (1/A^)(1 + ni • A) and 
P2 — (1/A^)(1 + ^2 ■ A) are orthogonal if 

where 6 is defined by fii ■ n2 — cos 9. Note that for N — 2 this reduces to the well-known fact that for two-state 
systems, the orthogonal states are represented by antipodal points on the Bloch sphere. 

The first condition in Eq. ® implies that the coherence vector must have unit magnitude. This restricts the set of 
vectors to those that lie on the surface of the unit sphere S'^^^. The second condition restricts the set of allowable 
rotations to a proper subset of the group SO{N^ ~ !)• The equations are non- linear and give a set of constraints which 
restrict to the manifold CP^^^ having 2N — 2 dimensions. The second condition is also related to the positivity of 
density operators, a fact which is discussed further below. 
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B. Mixed States in A'^-dimensions 

The mixed state density operator in A^-dimensions can be written in the same form as the pure state case: 

in{n- 1) 




n 



•A . (9) 



with n ■ n < 1. However, unhke the case for a two-state system, there are more constraints on the coherence vector 
for dimensions greater than two for the Hermitian matrix here to represent a positive, semi-definite operator. This 
wiU be given in Section IlIII 

C. Casimir Invariants 

The Casimir operators are invariant operators constructed from the Lie algebra elements. In particular, they form a 
maximal set of algebraically independent elements of the center of the algebra, formed by homogeneous polynomials in 
the generators. A very general discussion may be found in [l^, and were first constructed in (9|. General expressions 
for these are given in Appendix^ Here we note that the values of these operators can be determined by their relation 
to the trace invariants. For example, let us consider a density matrix, p. For all p 

Tr(p2)^l(l + (iV-l)n.n). (10) 

The quantity n-n is the value of the quadratic Casimir operator (see Appendix^, which we refer to as the quadratic 
Casimir invariant. An example of the quadratic Casimir operator is the total angular momentum operator. The 
Casimir invariants are unchanged by unitary transformations on the density operator. Similarly, 

Tr(p3) = -i-[l + 3(A^-l)n-n+(7V-l)(Ar-2)(n*n)-nl, (11) 

is clearly invariant under unitary operations. The quantity n-kft ■ ft is the cubic Casimir invariant. In the appendix 
we give the expressions for Tr(p"), n < 9. One may then recursively find higher order Casimir invariants and show 
that they are indeed unchanged by unitary transformations. The trace invariants, Tr(p"), here were discussed in [T^ 
where some discussion of the local unitary invariants were given for GHZ states. 

D. Constructing Local Invariants 

We can now construct a set of quantities which are invariant under local unitary transformations. These invariants, 
like the Casimir invariants are a distinguished set. Clearly local unitary operations preserve the Casimir invariants of 
the marginal density operators. However, in this section we discuss invariants associated with the correlation matrix. 

As an example, consider the quadratic Casimir invariant 

C2 — ft ■ n. (12) 
A two-qubit density operator can be expressed in a tensor product basis as 

p = ^(l® ]l-|-nA-CT«)H-]l(g)nB • (7 -f Cy-CTi «) CTj), (13) 

Note that local unitary transformations on systems A and B, denoted Ua and Ub, conserve ft^ ■ ft^ and fiB ■ fts 
respectively. This can be seen as follows. 



where R e 5*0(3). We can therefore rewrite 



UAn^^(J,u\^n^^^I^^aj, (14) 



vf^Ri^n^, (15) 
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and note that ua ■ nA = 'fnA ■ iriA since the transformation is orthogonal. We also know that the set of all unitary 
transformations acting on the composite system will be a subset of the matrices in S'0(15). This implies that 



nA-nA + UB • + CijCij , (16) 



is also a conserved quantity. However, we may want to ask what quantites associated with the correlation matrix, 
Cij, are conserved under local unitary transformations. The correlation matrix has rows and columns labeled by the 
indices i and j respectively. Now consider the vector formed from the elements in each. Examining Eq. (|14(l . we 
see that the magnitude of these vectors, is conserved by Ua- Similarly, the magnitude of the vectors formed by the 
columns is conserved. We may express these relations as, 

UACij cTj (g) (Jj = C'lj (Ji(E)aj, (17) 

where = RlCij , implies 

CijCy = ^ CyCy . (18) 

i I 

Similarly for Ub acting on the vectors formed from the columns of Cij . Therefore under local unitary transformations 
of the form Ua Ub, the following quantity is conserved, 



(19) 



More generally, we may determine conserved quantites formed from the correlation matrix which are analogues of 
the Casimir invariants. For the cubic Casimir invariant, for example, the following quantity is invariant under local 
unitary transformations, 

^ ^ dijkdlranCilC jmCkri' (^^) 
ijk^lran 

Similarly, we could construct invariants for systems of arbitrary dimension as well as systems with any number of 
subsystems. 

The number polynomial invariants under unitary transformations grows rather rapidly with the dimension of the 
system under consideration |l5j . One might suppose that only a subset is required for constucting entanglement 
measures given that, for example, the square of the concurrence |la,ll3 ^^r two qubits (see Section llV CI) is constructed 
from only three quantities which are invariant under all local unitary transformations. Here we have given a subset of 
local invariants which may well be useful for many quantum information processing tasks. The set of invariants given 
by Makhlin Ts] (see also 0,E9) to determine equivalence under local unitary operations is larger than the number 
of Casimir invariants, which are included as a subset, and are a complete set for determining the ability of two density 
operators to be transformed into one another by local unitary transformations. However, since the concurrence and 
I-concurrence poj do not rely on this large set of invariants, one may expect, generally, the number of invariants 
needed for the construction of entanglement measures may be far less than the number required for other purposes, 
such as local unitary equivalence. 

We have now shown that a density operator can be parameterized in terms of a set of traceless, orthogonal, Hermitian 
matrices and have constructed associated invariant quantities. Our next goal is to give positivity constraints for the 
density operators that determine the allowable sets of coherence vectors n. 

III. CHARACTERISTIC POLYNOMIAL/POSITIVITY 

In this section the characteristic polynomial of the density matrix is expressed in terms of the trace invariants and 
the Casimir invariants. 

A. The Characteristic Polynomial 



In this subsection we express the characteristic polynomial in several different ways in terms of invariants of the 
group. Consider an n x n complex matrix A of arbitrary dimension with eigenvalues pi. The characteristic equation 
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for the matrix can be written as (for a similar expression, see |21|) 

det{A - Al) = A" - 5iA"-i + S-sA""^ - +... + (-1)"5„ = 0, (21) 
where the Sk are the symmetric functions given by |22| 

fc 

Sk^ E n^''.- (22) 

i<n<---ifc<A'j=i 

These can be written in terms of [Tr(p")]'" as 

^i=Tr(A), ^2 = (l/2)[Tr(A)5i-Tr(A2)], (23) 

and 

Sk = {l/k)[TY(A)Sk-i - TT{A^)Sk-2 + ... 
+ (-l)"-iTr(^")5,_„ + ... 

+ {-1)''-^Tt{A''-^)Si + (-l)'^-iTr(A'=)]. (24) 

This can be proved using the fact that 



N 



--J2iN;m,,m,,...,mM)prPr ■■■P'm'^ (25) 



where {nik} is a set of integers such that X^feli ™fe ~ ^^'^ 



(iV;mi,m2, ...,mM) = ; ; r. (26) 



B. Positivity 



For a given set of real numbers {ni,n2, ...,nAr} € M^, we would like to know when the set will represent a valid 
density operator of the form Eq. ijU)). It is clear that the right hand side of Eq. © has trace one and is Hermitian. 
However, the positive semi-definite property is less trivial. 

Theorem: For a Hermitian matrix p = (1/7V)(1-|- y^(]V(]v"^~T)72)~n . A) to represent a positive semi-definite operator 
it is necessary and sufficient for Sk > for all k. 

Sketch of proof: Since the matrix p is Hermitian, all eigenvalues of the operator are real. This implies that the 
coefficients of the characteristic polynomial are real. They are also non-negative if and only if the signs of the 
coefficients of the characteristic polynomial alternate. In fact, the number of positive roots of the characteristic 
polynomial is the number of sign changes in the sequence of coefhcients (pages 124-5, 23J). □ 



1. Constraints on the Coherence Vector 



The set of inequalities 5*^ > characterizes the region of permissible vectors which represent valid, i.e., positive 
semi-definite, density operators. The first few of these conditions, given directly in terms of the coherence vector, are 
as follows. For a normalized p, 



5i =Tr(p) = 1. 



(27) 
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Here we adhere to the conventions set forth in Sections III Al and III Bl Using the symmetric parts of the traces, denote 
Tigyjn given in Appendix IBl 



^2 
^3 

Si 



l[(Tr(p)f-(Tr(p2))]^^[l-n.n], 



1 {N- 1){N~2) 



6 

1 (iV-l)(iV-2)(7V-3) 
24 

1 — 6n • n + 8n ★ n • n + 



[1 — 3r? • n + 2n * n • n] , 



(28) 
(29) 



6(iV- 2) 



-n-k n ■ n-k n 



(30) 



Higher order invariants can be calculated using the material from the Appendices in a straightforward albeit 
somewhat tedious manner. Note that if the two requirements for a density operator to be a pure state are met, 
n ■ n — 1 and n-kfi — n, then 5*2 through S4 (as well as all higher 5*^) vanish, indicating a characteristic polynomial 
with the solution, one non-zero eigenvalue. The trace being one then demands that this eigenvalue be one. 

It is also noteworthy that two density operators have the same Casimir invariants if and only if they have the 
same eigenvalues. This follows from the fact that two density operators have the same Casimir invariants if and only 
if they satisfy the same characteristic equation. An entanglement measure based upon an entanglement monotone 
for a bipartite pure state must be a function only of the eigenvalues of the marginal density operators T^. This 
relation between Casimir invariants and eigenvalues implies that any entanglement measure based on an entanglement 
monotone may also be expressed as a function of the Sk or Casimir invariants of the marginal density operator. 



C. Symmetric Functions and Casimir Invariants 



The quantities appearing in the Sk are combinations of the Casimir invariants. This relationship is noteworthy 
for reasons other than those just stated. Casimir invariants can be used to determine degeneracies in the orbits and 
emphasizes the relation to the physical system and Casimirs invariants are conserved quantities used as labels for 
quantum states. To illustrate the ability of the Casimir invariants to provide information about the degeneracy of 
the spectrum, we will use the three-state system as an explicit example and then give a brief discussion of four-state 
systems. 



1. Casimir Invariants for a System with Three States 

Since the eigenvalues are invariant under unitary transformations, we can discuss the interpretation of the Casimir 
invariants in terms of a diagonalized density operator. In three dimensions a common basis for the traceless, diagonal 
3x3 Hermitian matrices are the Gell-Mann matrices |2^. In this basis, we denote the two linearly independent, 
traceless diagonal matrices as 

^ /I \ 1 /I o\ 

As - -= 1 , A3 = ^ -1 . 

For a mixed state, we may write the diagonalized form as 

f ai \ 
Pd = as , 
V 03 / 
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where '^^cii = 1- Expanding this using 



Pi 



1 + ^ f V3 A3 + As 



92 = ^ 



P3 



yields a density operator of the form 



1 







1 







1 



(31) 
(32) 

(33) 



1 + V3((ai V3/2 - a2%/3/2)A3 
+ (ai/2 + a2/2-a3)A8^ 



(34) 



The coherence vector is given by 



= (0, 0, aiV3/2 ~ a2\/3/2, 0, 0, 0, 0, ai/2 + 03/2 - 



03) 



Since this is a positive semi-definite, Hcrmitian matrix, the density operator formed by pm = UpdU^ ^ ^(^'^ V^Un- 
XU^) is also a positive semi-definite, Hermitian operator. With the appropriate restrictions on the coefficients, we 
may par ameterize all three-state density matrices (and a direct generalization for higher dimensional systems) in this 
way |25ll26l|. 

For three-state systems, the following two quantities are two independent Casimir invariants which, in terms of the 
coherence vector, are given by 



72 ■ n — C2, 



n-k n ■ n — c^, 



(35) 



The first is the quadratic Casimir invariant of the group and the second is the cubic Casimir invariant of the group 
(see also |1,0|). The generic orbits are given by 21], 



fi ■ fi — C2, and ft-k fi ■ n = ^ C2. 

The values of C2 and C3 are unchanged, i.e. invariant, under unitary transformations of the density operator. 
The square of the coherence vector is 



(36) 



We may also calculate 



fi ■ fi — a\ + a\ + — 0102 — aia^ — 0203 < 1. 



n-kfi-n — 0^+02+03 + 6010203 

— (3/2)(oia2 + O2O1 + 0^03 + 0303 
+0301 + 0302). 



Note that —\nr < n-kn ■ n < \nr since 



27, 



{n-kfi-n)^ - \n\^ = ^(oi - 02)^(01 - 03)^(02 - 03)^ > 0. 



Degenerate eigenvalues imply the following relations. 
1. if oi = 02 



n • n = (oi — 03)^ and n • n = — (oi — 03)'^ 



(37) 



2. if 02 = 03 



n ■ n = (oi — 03)^ and n-kn - n — (oi — 03)'^. 



(38) 
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3. if fli = 03 

n ■ n = {a2 ^ 0,3)'^ and n -k n ■ n ^ {a2 — a^)^ . (39) 

Therefore, when the two eigenvalues are degenerate, n-kn- n cx When the two degenerate eigenvalues are greater 
than the third, the quantity nkn- n is negative and when they are smaller, ft-kn- ft is positive. Thus by investigating 
the values of the Casimir invariants, we are able to extract information about degeneracies in the spectrum. These 
degeneracies correspond to invariant subspaces since an eigenvalue subspace spanned by degeneracies is invariant 
under unitary transformations on that subspace |25| . We next comment briefly on the four-state and general cases of 
identifying degeneracies. 

2. Higher Dimensions 

For iV-state systems, there are TV— 1 Casimir invariants. This is the rank of the group of transformations, SU (N) on 
the space of density operators, and corresponds to the number of elements in a complete set of commuting operators. 
Each Sk, when expressed in terms of the coherence vector, will contain a term of the form {n-k)^'-^~^^n ■ n which is 
absent from Sj, j < k. In the previous section it was shown that a degeneracy in the spectrum of the density operator 
was manifest in the values of the Casimir invariants. When a degeneracy in the spectrum exists, an added symmetry 
of the density operator under a subgroup of the group of all unitary transformations exists. This will determine a 
relation between the Casimir invariants, and thus reduce the number of independent polynomial invariants. 

Let us discuss the example of four-state systems. If the density operator for a four-state system has the following 
spectrum, (a,b,b,b) then the each of the four Casimir invariants are proportional to powers of \n\ (ci cx If 
the spectrum is {a,a,b,b), then all Casimirs are zero except the quadratic. Spectra of the form (a, 6,c, c), or non- 
degenerate spectra are not as easily idenitfied by their Casimir invariants. However, there exists a readily available 
program, Macaulay, which can check the independence of the invariants, thereby determining the degeneracies. Of 
course, if the spectrum is completely degenerate, then all Casimirs vanish since n = for the completely degenerate 
case. The advantage of obtaining this information through the use of invariants is that one may not always solve 
directly for the eigenvalues of a matrix, but the Casimir invariants may still be obtained. 

For the convenience of the reader, the Casimir invariants are given in terms of the Lie algebra elements in Ap- 
pendix^ In Appendix IbI we give the trace formulas from which these can be calculated and the coefficients of the 
characteristic polynomial can be found. 

Note that a map from a density operator to a density operator may be expressed as an afhne map, 

n^n' =Tn + t, (40) 

where T is a matrix and t is a translation. The positivity of the mapping is determined by the positivity of the density 
operator formed by n' ^ . 

IV. EXAMPLES 

In this section we give the following results. First, we show how the positivity of the Sk restrict the coherence 
vector for two particularly interesting examples, nkft - n — —\n\^ and inversion. This gives, in terms of the coherence 
vector, the same bound obtained by Rungta, et al. po| on the ability to construct a "universal inverter." Second, we 
show that the positivity of the density operator of two qubits can be determined by the positivity of 5*3 and S4 for 
the general case and for the Werner state. Third, we present an alternative derivation of the three-tangle of Coffman, 
Kundu and Wootters using the coherence vector description. 

A. Inversion of the Coherence Vector 

Here we show that, due to positivity requirements, the limit n-kn - n = — \nf' cannot be reached for certain n. This 
follows from the positivity requirements Sk>Q and restricts the set of positive maps for the set of density matrices. 
An example of this is the universal inverter and universal NOT gate. 
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1. Universal Inversion 



The universal inverter and universal NOT gate j20j are related to a mapping of the form 

p^l-p, (41) 
which is positive but not completely positive. In terms of the coherence vector representation, 

P - ^ (l(iV -D-cn-X)^ {t - -J-n . a) , (42) 

where c = ^JN{N — l)/2. Thus, up to an overall constant, the mapping corresponds to a change in sign of the 
coherence vector and a reduction of the magnitude of the coherence vector. 



2. Inverting the Coherence Vector 



We might ask if there exists a physical map which will properly invert the coherence vector. (Inversion of the 
coherence vector as a possible generalization of the concurrence |1^ iHiJ was studied by Rungta, et al. i2Qj.) This 
would be of the form 



P= -^(1 + cn- A) -> p = 
However, this is not positive. To see this, consider the matrix 



1(1 -en- A). 



(43) 



1 



/ a 
a 

\0 ... 







-{N-l)aJ 



(44) 



For this matrix fi ■ n = a? and ni^ n ■ ft — a^. This gives the symmetric polynomial 

S3 cx 1 — 3n ■ n + 2n-k ri ■ n = 1 — 3a^ + 2a^ . 
This function of n is minimum when n-kn ■ n = — InP < so that 



5*3 cx 1 — 3n-n — 2|n*n-ri| = 1 — 3a^ 



2a^ 



(45) 



(46) 



For this to be positive, a > 1/2 showing that for certain n the limit n -k n ■ n — —\n\^ cannot be obtained. This is 
unlike the case of a Hamiltonian, or general Hermitian matrix, where it is acceptable to have n-kfi ■ ft ~ — \fi\'^. For a 
system with three states, and no zero eigenvalues, S3 is the non-zero determinant of the matrix. 
For higher dimensional systems the requirement that p in Eq. (|44|l be positive corresponds to 



1 



> a > -1. 



N - 1 

Now if we ask for an inversion map which is positive, we seek a mapping of the form 

1 



P 



N 



{hi- eft- X). 



Choosing an operator of the form Eq. (|44|l . for the map to be positive we require 

b> {l-N)a> {N -1). 



(47) 



(48) 



(49) 



This is the condition found by Rungta, et al. 20] for positivity and restricts inversion to a map of the form in Eq. H41|l . 
This is a condition on the positivity of the determinant which we have shown is Sn for an iV-state system. 
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3. Three-State Example 



For example, let us consider a three-state density matrix of the form 



0.15278 0.036084- i0.06250 -0.072169 + i0.12500 

0.036084 + i0.06250 0.23611 -0.25 | . (50) 

-0.072168-^0.12500 -0.25 0.61111 

Using jii = (\/3 /2)Tr(pAi), direct calculation gives 

5*3 cx 1 - 3(0.666)2 ^ 2(0.666)^ 
However, when 0.666 — > —0.666 then 6*3 < showing that inversion is not a positive map for this density operator. 



B. Two Qubit Entanglement 



In the next subsection IjIV B 2|) the example of the Werner states for two qubits is investigated. This mixture of 
a completely mixed and singlet state is separable if and only if the partially transposed density operator is positive 
semidefinite according to the Peres-Horodecki criterion In this case, S'3 and 6*4 determine positivity. This will 

be shown using the coherence vector representation. 



1. A Basis for Two Qubits 
Let a basis for the Lie algebra of SU{4:) be given by 

{K}lio = Wz<^'yj}lj=o, (51) 
where Aq = 14 and ctq = ^2- The labels correspond in the following way, 

A,, i = 0,1,2,3 -^ai(g)l, i = 0,1,2,3, 

Xi, i = 4, 5,6 ^ ^ 10(7^, z = 1,2, 3, 

A,;, i = 7,8,9 ^ ^ cTi ®cri, i = 1,2,3, (52) 

Ai, i = 10,ll,12 ^ -ij 0-2 «)(Ti, z = 1,2,3, 

Ai, « = 13, 14, 15 ^ -ij 0-3 (Ti, z = 1,2,3. 

This forms an orthogonal basis with respect to the trace and has normalization given by 

Tr(A,Aj) = 2%. (53) 
The non-zero, totally symmetric d-tensor components in this basis are given by: 

—F= — dl,4.7 = dl.5,S — di.eg — (i2,4.10 = ^2,5.11 — ^2,6.12 

— 1^3,4,13 = C?3,5,14 = '^3,6,15 = — C?7,ll,15 = — C?8,12,13 

— 1^7,12, 14 — — ^9,10,14 — C^S^O.IS = 1^9, 11, 13- (54) 



2. Werner States: A Case Study 



Under partial transpose of the first subsystem in the density operator, only elements rt2, nio, nn, ni2 change sign (in 
the given basis Subsection IIV B l]) . Therefore under the partial transpose, one may readily determine which elements 
of the products n* n ■ n and n* n ■ n-k n change sign. 

The inequalities S'3 > and ^4 > depend only on the non-local invariants of the system since S2 does not change 
and the local invariants which have the same form of S2 also do not change. This shows that the negativity arises 
in the nonlocal invariants (as they should). As noted before, the partial transpose is positive since it preserves local 
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positivity, but is not completely positive. Although this is a low-dimensional example and the higher order Sk become 
more complicated as the k increases, such an analysis might lead to ways (e.g. numerical and/or analytic searches) 
for identifying positive, but not completely positive maps which may witness entanglement. 

To clarify the discussions above concerning the positivity of the coefficients of the characteristic polynomial, we 
give an example of the calculation for the Werner state of two qubits. The Werner state for two qubits is given by 



1-x 



Pw = 



1 + xS, 



(55) 



where < x < 1 is real and S is the singlet state 




1-10 
0-110 




(56) 



Therefore when a; = the state is separable and when x = 1 the state is maximally entangled. We may rewrite this 
as 



1 X 
Pw = ^1 - ^(o-x ® Ox 



1 


4 















l+x 

4 


X 

2 










X 

2 


l+x 
4 





\ 











\ — x 
4 / 



(57) 



The partial transpose condition (Peres-Horodeckis) '5, 's'l is equivalent (up to a local unitary transformation) to the 
inversion of the coherence vector, which is also known as spin flip or inversion. In terms of the coherence vector for 
the combined system, if we write the density operator in terms of the basis given in the previous section. 



Pw — Pab 



1 



1 + Ve n • A ) , 



(58) 



the partial transpose corresponds to n2 — + — ??2, "^lo ~* "'^lOj '^ii ~* —"^11,^12 —ni2- Calculating the coefficients 
of the characteristic polynomial, we find Si{pab) and S2{pab) are unchanged under this transformation. However, 



AB) 



^1(1 



_^ 2x3) 
(1 - 3x2 _ 22,3)^ 



(59) 



and 



S^ipAB) 



^ ) (1 - 6x2 + g^3 _ 3^4-) 



(1 - 6x2 _ ^^4y 



(60) 



This partial transpose condition implies that the density operator is separable if and only if the partially transposed 
density operator (or the spin flipped density operator) is positive semi-definite. Here we see that the coefficients have 
following possibihties for sign changes. For 1/3 < x < 1/2, 5*4 < 0, 5*3 > 0, and for x > 1/2, £'4 < and 5*3 < 0. 
However, in each case there is only one change in sign for an Sk and therefore one negative eigenvalue. 



C. Distributed Entanglement 

Coffman, Kundu and Wootters [23| have studied "distributed entanglement" which concerns the entanglement of 
various subsystems of a tripartite qubit system. One of their main results is the description of entanglement of a 
pure state of three qubits which is not expressible in terms of two-qubit relations. Here we wish to streamline their 
argument using the material presented above and thus derive by alternative means the "tangle" of three qubits. 
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Consider a pure state of three qubits for systems we label A, B, C. We will write the density operator in a tensor 
product basis, 

PABC = ^(l^il^l + nA-CT^l^l+l^ns-CT®! 
8 

+nAc • 1 (8i (? + ubc • 1 (? (? 

+nABC ■ ^ ^ (61) 

where uab ■ (J ® a = {nAB)ij<^i ^ dj etc. 

Since pabc represents a pure state, the marginal density matrices, e.g., pab — Trc(/3ABc) has only two non-zero 
eigenvalues, so that the square of the concurrence may be used to write 

Cab ~ (-^1 ~ -^2)^ = + A2 — 2A1A2 

= Tt{pabpab) - 2A1A2 < Tt{pabpab), (62) 

where Ai and A2 are the square roots of the eigenvalues of pab Pab- The matrix pab is defined by pab — ® 
^vPab^^v ®^y- 

At this point our argument will differ from that of ,27:]. Since this is a pure state, the Schmidt decomposition can 
be used to choose a preferred basis for subsystems AB and C . The reduced density matrices may be rewritten as 
(using an unnormalized coherence vector) 

PAB = Trc(/9ABc) = ^(1 + ^AB ' A), (63) 

where rfiAB = {nA^nBTnAB) and 

pc = T^^AB{pABc) = \{^ + nc-^)- (64) 

Then, by the Schmidt decomposition these two have the same eigenvalues. Therefore they satisfy the same charac- 
teristic equation which will have only one non-trivial Sk [Si — Tr(p) — 1), that being 

Siipc) = S2{pab), (65) 

which implies 

1 1 

-{1 + mAB ■ mAB) = -^{'i^ + nc ■ nc). (66) 

Therefore 

nAB ■ nAB = 1 + 2?ic • nc - ua ■ ua - ub ■ ns (67) 

Noting that 

Tr(pAsPAs) ^^{l-UA-nA-nB ■nB + uab ■ nAB), (68) 

^i:{pabPab) = 2 (1 ~ ftA ■ ftA - ub ■ fiB + fic ■ nc). (69) 
This is completely equivalent to the results in Eqs.(7) and (8) of [23|, the latter is repeated here: 

Tt{pabPab) = 2(detpA + detps - detpc)- (70) 
This is needed to derive the "first main result" of 27]: 

Cii3+Cic<4dctp^, (71) 

where we have used Eq. (|62|) . 
At this point, we can calculate 



we can use Eq. (|67|) . to write 



4\/ S2{pabPab) — tabc = CIa)bc ~ ^\b ~ C^c- 

This quantity describes the three-way entanglement of the three qubits and was shown in p7j to be invariant under 
the permutation of the qubits. 
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V. CONCLUSION 



The identification of positive but not complete positive maps fias recently become an active area of research due 
to the restrictions it places on physically realizable quantum transformations and the question of entanglement 
of quantum systems 0, 01 ■ To aid in the study of such transformations this paper has presented a representation 
of the density operator in terms of traceless, Hermitian, orthogonal matrices. We then showed that the Casimir 
invariants of generalized coherence vector for density operator could be calculated directly and information about 
degeneracies in the spectrum of the operator could be obtained for some particular cases. It should be noted that 
we have given a representation of the density operator in bases, but the expressions of the Casimir invariants and 
symmetric functions do not depend on the choice of the set of traceless, Hermitian, orthogonal matrices in the basis. 
The region of positive semi-definite density operators is determined by the necessary and sufficient conditions, Sk > 0. 
The Sk were expressed in terms of the coherence vector and Casimir invariants. The positivity conditions given here 
not only indicate whether a density operator has all positive eigenvalues, but it also indicates the number of positive 
eigenvalues in terms of the number of sign changes of the sequence of coefficients Sk- 

Superoperators which map Hermitian operators to Hermitian operators will preserve the reality of the eigenvalues. 
Since the eigenvalues are real, the coefficients of the characteristic polynomial must alternate in sign if the eigenvalues 
are to be positive. Therefore changes in the signs of the Sk can indicate positivity or non-positivity of maps of the 
density operator. Given the expressions in this paper, this statement may be utilized directly given an affinc map of 
the coherence vector. 

It is interesting to note that the "measure of purity" of a density operator has arisen in several contexts. Consider 
a pure state, bipartite density operator. The generalized concurrence in poj is simply related to the purity of the 
marginal density operator. If pA is the marginal density operator, then the concurrence is proportional to S2{pa) 
which is a measure of the purity of the density operator. The state pA is pure if and only if S2{pa) is zero. The state 
is "less pure" if this quantity is larger. This measure of purity is also used in the optimal decompositions discussed 
in p^ . One might consider generalizations of the "measure of purity." Certainly if Si (equal to one when the matrix 
has unit trace) and S2 are the only non-zero coefficients of the characteristic polynomial, then ^2 is a "good" measure 
of purity. However, if ^2 and S3 are both non-zero, then the purity should be measured by two quantities since pure 
states necessarily have both quantities equal to zero. States that are closer to being pure are those with smaller values 
of these two quantities. Similar arguments can be made for the higher dimensional Sk- One might then consider a 
generalization of measures of entanglement which rely on this modified set of "measures of purity." 

The set of algebraic equations given by Sk > give a set of geometric constraints on the spaces of allowable 
coherence vectors. This may motivate further exploration of techniques from algebraic geometry which has already 
been found useful by Miyake |2^J for describing pure state separability. 

Due to the generality of the arguments here and the connections made between Casimir invariants, algebraic 
geometry and positivity, we believe this work provides useful relations and insights into the structure of positive 
operators. We also hope that it will aid in identifying positive, but not completely positive maps. 



APPENDIX A: CASIMIR INVARIANTS 

Here we give expressions for the Casimir invariants of a Lie group. For a discussion see |13|. 

The Killing form Gab gives the metric gab on the vector space. This will determine the quadratic Casimir invariant 

JV 

C2^Y1 dabX^^', (Al) 
a, 6=1 

where N is the dimension of the vector space (N = — 1 for SU{n) groups), and A € 'C(G). Note that gab oc 
d f'^'^df^'^c is invariant, symmetric tensor. To find other invariant, symmetric tensors, one forms 

TT{adxai o ad\a.2 o • • • o ad\an ) = 

N 

E ■ ■ ■ r"-''"\r"\ (A2) 

f)l,f)2,---,f)n = l 

One can express the Cubic Casimir invariant in terms of the totally symmetric tensor dabc, 

N 

C3= J2 dabcX^'X^'X'- (A3) 

a,b,c—l 
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Generally these higher order invariants can be expressed in terms of the symmetric tensor as 

Cm ~ ^ ^ ^016162 

ai,a2,---,0^m-3 
61 ,b2,---,bm 

X'^aia263'^a2<l3b4 • ■ • '^0^-2 0^-3 5m-2^am-36Tn-li>m 

xA^'^A''^ •••A'"'" (A4) 
We list the first few here in order to be explicit and to enable the development of the pattern. 

C4= J2 da,b,b2da,b,b,X'"X'''X'"X'" (A5) 

0,1 ,61 ,62 ,63 .64 

C5 = ^ <ioibi62t^aia263C^a2a465A''^A''^ • • • A^^ (A6) 

61, 62, 63, 64,65 

C'e — ^ ^ '^aibi 62^^010263 f^a2a364'^a36566 

01, a 2, "3 
61,62, ...,66 

xA^'^A''^ •••A''^ (A7) 

Of course the ones that are immediately interesting are C2, C3, C4, Cg, Cg for the purposes of embedding 2 qubits 
into a 4-state system, a 2-state and 3-state system into a 6-state system and the embedding of a two 3-state systems 
into a 9-statc system. These are useful for examining quantum control for two-qubits and entanglement issues for a 
two-qubits, a qubit and a qutrit, and two qutrits. 

The above relations can be expressed in terms of adjoint vectors and particular products. We introduce this 
notation here since it has its own manipulation rules that make it easier to calculate quantities of interest. Note also 
that since the fabc and dabc tensors are obtained by taking traces of products of elements with anticommutators and 
commutators respectively, they are easily calculated by analytic methods on a symbolic manipulation program such 
as MATHEMATICA. These relations are 

/„6c = Tr([A„,Afe]Ae), 

and 

dabc = Tr {{Xa, Xb}Xc) . 

The difference between upper and lower indices is not important if we are considering SU{n). 



APPENDIX B: TRACE FORMULAS 

1. Symmetric Traces of Basis Elements 
Here the first few examples of the trace formulas have been given. 

Tr(AiA,) = 2% (Bl) 

Trsj^m(AiAjAfe) = 2dijk (B2) 
4 

Trsj,TO(AiAjAfcA;) = —SijSki +2dijmdmki (B3) 
4 

Trsj/m,(AiAj AfcA;Ag) = —{Sijdklq + 6kldijq) + 2dijmdklndmnq (-^4) 



Trsi/m(AjAjAfeA;Aq,As) — -;^5ijSkiS, 



23 
4 

(B5) 
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23 



~^ j^i^qsdijmdklndmnu ~l~ ^ijdklmdqsndmnu ~l~ ^kldijmdqsndmnu) 
+^dijmdklndqstdmnt (B6) 



2^ 

H~ yy-2 (.^ij^kldqstduwt ~l~ ^ij ^qs^klnduwri H" ^ij^uvjdqsndhln 
+ SklSqsd 

ijnduwn H" ^kl^uwdijndqsn H~ ^uw^qsdklndijn) 
i^^ij dklmdqstduwvdtvm H" ^kldijmdqstduwvdtvm 

-\- Sqsdijrndklndrrvnvduwv ~\~ ^uwdijmdklndmnvdqsv) 
~\~ 2 ^ J rndkln dm np dqst du w v dtvp 

(B7) 

2^ 

'^■'^sym (•^i'^j •^fc'^^ '^Q'^S'^U'^W'^y ) — yy*3 i^ij^kl^qsduyw H" ^ij^kl^uwdqsy ~\~ ^ij^qs^uwdkly ~\~ ^kl^qs^uwdijy^ 

23 

~l~yy-2 {^ij^kldqstduvwdtvy H" ^qs^uvjdijmdklndmny 
~\~ ^ij^qsdklnduvvjdnvy H" ^ij^uwdklndqstdnty^ 

23 

H~yy-2 (.^ijdklydqstduwt ~l~ ^kldijydqstduwt 

~\~ ^qsdijmdklmduwy H" ^uvjdijmdklmdqsy^ 

22 

~t~ ^qsdijmdklndmnpduwvdpvy ~\~ ^uwdijmdklndmnpdqstdpty) 
~\~^dijijf.dklndmnpdqstduwvdtvxdpxy (■^^) 



~l~ J- (y^ijdklndqstduvwdtvxdnxy H" ^kldijmdqstduvwdtvxdnixy 



2, Symmetric Traces for the Density Operator 

For the density operator these translate to (again only the first four are given): 

'IV(p2) = l[i+(7V-l)ri.n] (B9) 

Tr(p3) = _L[i + 3(7V-l)fJ.fi+(7V-l)(7V-2)(n*n)-n] (BIO) 

Tr(p4) = i3[l + 6(Ar-l)fJ.rJ 

+4(7V - 1){N -2){n-kn)-n+{N - if {ft ■ nf 

+ {N-l){N-2f{ni.n)-{n*n)] (Bll) 

Tr(p5) = -^[l + 10(7V-l)n-n+10(7V-l)(7V-2)n*n-n 

+ 5{N - l)2(fj . fj)2 + 5(iV - 1)(A^ - 2)2(ri • (fiTt n) 
+ 2(7V- l)2(Ar-2)(n-n)(n*n-n) 

+ {N -1){N -2fn*n*n*n-n] (B12) 
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Tr(p^) = -^[1 + 15{N -l)n-n + 20{N -1){N -2)ni<n-n 

+ 15{N-lf{n-nf + 15{N-l){N-2f{ni.n) -{n-kn) 

+ 12(iV- lf{N-2){n-n){ni.n-n) 
+ 6(iV -1){N -2f{n-kn-kn-kn-n) 

+ {N-lf{n- nf + 3(iV -lf{N-2f{n-n){n-kn-nirn) 

+ {N -1){N -2f{n-kn-knf] (B13) 

Tr(p'^) = -^[l + 21(iV - l)n • n + 35(Af - l)(Ar - 2)n*n • n 

+ 35(iV- l)2(n.n)2 + 35(Ar_ i)(7V-2)2(n*n) • (n*n) 
+ 42(Af- l)2(iV-2)(n-n)(n*n-n) 
+ 21(7^- 1){N -2f{n-kn-kn-kn-n) 

+ 7{N-lf{n-nf + 21{N-lf{N-2f{n-n){n*n-n*n) 

+ 7iN -1){N- 2)^{n irfi^nf 

+ 3{N -lf{N-2){n- nf {n-^n-n) 

+ 3{N-lf{N-2f{ni^ni^n*n-n){n-n) 

+ {N - 1){N - 2f{nikn - n){ni<n- n-kfi) 

+ {N -1){N -2f{n-kn-kn-kn-kn-kn-n)\ (B14) 



Tr(p^) = -^[l + 28(7V- l)n-n + 56(7V- l)(7V-2)n*n-n 

+ 70(7V - lf{n ■ nf + 70{N - 1){N - 2f{n-kn) ■ {n*n) 
+ 112{N-lf{N-2){n-n){n-kn-n) 
+ 56(iV - 1)(7V - 2)^(n ★ n ★ n ★ n • n) 

+ 28(A^ - lf{n ■ nf + 84(iV - lf{N - 2f{n ■ n){n*n ■ n-kn) 

+ 28(7V - 1)(7V - 2)*(n ★ n ★ nf 

+ 24{N-lf{N-2){n-nf{n*n-n) 

+ 24(Ar- l)2(iV-2)3(H*n*n*n-n)(n-n) 

+ 8{N - 1){N - 2f{nkn- n){nkn ■ ni^n) 

+ 8{N - 1){N - 2f{n-kn-k n-kn -k n-kn ■ n) 

+ {N- lf{n ■ nf + 6{N - lf{N - 2f{n ■ nf {n-kft ■ n-kn) 

+ A{N - 1)'^{N -2)'^{n-n){n-kni<ni<ni<n-n) 

■j- {N — 2)^{n-kn-kn-kn-kn-kn-kn ■ n)] (B15) 
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Tr(p^) = — [l + 36(iV-l)n-n + 84(7V-l)(iV-2)n*n-n 

+ 126(iV- l)2(n- + 126(iV- l)(iV - 2)2(r7;*n) • 
+ 252(iV - l f{N - 2)(r7; • n){n*n ■ n) 
+ 126(iV- 1){N - 2f{n*n*n*n - n) 

+ 84:{N-lf{n-nf + 252{N- if {N - 2f {n ■ n){n n ■ n -k ft) 

+ M[N -\){N -2f{nirnirnf 

+ 108(iV - lf{N -2){n- nf{n-kn ■ n) 

+ 108(iV - lf{N - 2)^(n ★ n ★ n ★ n • n){n ■ n) 

+ 36{N - 1){N - 2f{n*n- n){nirn- n*n) 

+ 36(7V - 1){N - 2f{n-kn-kn-kn-kni^n - n) 

+ 9{N - l)*{n ■ nf + 54(7V - l f{N - 2f{n ■ nf{n-kn- n*n) 

+ 36{N-lf{N-2f{n-n){n*n-kn*n*n-n) 

+ 9(-/V — 2)^(n*n*n*n*n*n*n-n) 

+ 4(7V - lf{N - 2)(n • nf{n*n ■ n) 

+ 6(7V- lf{N- 2f{n-nf{n-kn-kn*n-n) 

+ 4(iV - l)^(iV - 2f{n ■ n){n*fi- fi){n*n- n*n) 

+ 2{N - if {N - 2f (n ir n ■ n ir n){n ir n -k n -k fi ■ n) 

+ 4(7V - if {N - if {n ■ n){n k ni< n-k n ■ n-k fi-k fi) 

+ {N — 1){N — 2f {n * n n n i< fi -k n -k n -k n ■ n)] (B16) 
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